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Abstract. The correlation term, introduced in [T3] to describe the interaction between very 
far apart vortices, governs symmetry-breaking for the Ginzburg-Landau equation in R2 or 
bounded domains. It is a homogeneous function of degree (—2), and then for ^—symmetric 
vortex configurations can be expressed in terms of the so-called correlation coefficient. Ovchin- 
nikov and Sigal [13] have computed it in few cases and conjectured its value to be an integer 
I multiple of j. We will disprove this conjecture by showing that the correlation coefficient 

T-H ■ always vanishes, and will discuss some of its consequences. 
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^ ; 1. Introduction 

p ^ ■ The Ginzburg-Landau theory is a very popular model in super-conductivity [6]. Stationary 

^ . states are described by complex-valued solutions u of the planar equation 

\ where A; > is the Ginzburg-Landau parameter. The condensate wave function u describes 
, the superconductive regime in the sample by simply interpreting as the density of Cooper 

electrons pairs. The zeroes of m, where the normal state is restored, are called vortices. The 
^ . parameter k depends on the physical properties of the material and distinguishes between Type 
lO I I superconductors /c < ^ (in this normalization of constants) and Type II superconductors 

• \ In the entire plane M the parameter k does not play any role, as we can reduce to the case 
^ ■ A; = 1 by simply changing u into Supplemented by the correct asymptotic behavior at 

CN ■ infinity, the Ginzburg-Landau equation now reads as 

^' -^U = U{l-\U\'^) inM2 

^ ^ , |f/| 1 as l^l oo. ^ ' ' 
rS" 

\ The condition |?7| — 1 as |x| — oo allows to define the (topological) degree deg f/ of [/ as the 



winding number of f/ at oo: 



deg t/ = 7^ / d{argU), 

J\x\=R 



'\x\=R 

where i? > is chosen large so that \U\ > \ in \ 5/^(0). Given n G Z, the only known 
solution of (11. ip with deg U = n is the "radially symmetric" one Un{x) = '5'„(|a;|)(|^)" (in 
complex notations with x G C), where Sn is the solution of the following ODE: 

Sn + lSn-^Sn + Sn{l-Sl)=0 in (0, +Oo) 

S'„(0) = 0, lim S'„ = l. 
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Existence and uniqueness of Sn is shown in [HIT]. Moreover, the solution Un is stable for \n\ < 1 
and unstable for \n\ > 1 |Tl]. When n = ±1, the solution U±i is unique, modulo translations 
and rotations, in the class of functions U with deg U = ±1 and J^2{\U\'^ — l^dx < +oo [TU] . 

One of the open problems (Problem 1)- that Brezis-Merle- Riviere raise out in [3]- concerns the 
existence of solutions U of fll.ip with deg U = n, \n\ > 1, which are not "radially sjTiimetric" 
around any point. So far there is no rigorous answer, but a strategy to find them has been 
proposed in [12]. Formally, a solution U of (11.11) is a critical point of the functional 



Since £{^) = +oo for any C^— map \& so that l^'l — 1 as |x| — +oo and deg (\E') 7^ 0, 
Ovchinnikov and Sigal [TT] have proposed to correct S into 

where x is a smooth cut-off function with x = when \x\ < R and x = 1 when \x\ > 
R + R~^, and R >> 1 is given. Given a vortex configuration (a, n) = (^i, . . . , ax, • • • , nx), 
a C^— map so that — 1 as |a;| — > +00 has vortex configuration {a,n) if ai, . . . ,aK are 
the only zeroes of with local indices ni, . . . , nx, denoted for short as conf = {a,n)- Given 
Uq, Ovchinnikov and Sigal [12] introduce the "intervortex energy" E given by 

E{a) = inf{£ren(^) : conf $ = (a,no)}, 

and conjecture that is a critical point of E if and only if there is a minimizer U for i?(ao), 
yielding to a solution of (II. ip with conf U = (oojIIo) which is not "radially symmetric" around 
any point by construction. Letting da = niin \ai — aj\, the following asymptotic expression is 

estabhshed [T2] : 

K 

E{a) = J2^ren{Un,) + H{^) + Rem (1.2) 
i=i 

with Rem = 0{d~^) as da +00, where i/(a) = — tt ^ rijUj In |aj — aj\ is the energy of the 
vortex pairs interactions. When VH{a) = 0, the estimate in (II. 2p improves up to Rem = 

o{d-J). 

When VH{a) = (a so-called forceless vortex configuration), by choosing refined test functions 
the asymptotic expression (11.21) is improved [13] in the form of the following upper bound: 



K 

E{a) < ^ren(t/nj + H{^) - A{a) + Rem (1.3) 
i=i 

with Rem = 0{d~'^ + i?^^) as da — ?■ +00, where the correlation term A{a) is a homogeneous 
function of degree (—2) given as 

K K 

i=i i=i 

with V5j(x) = nj9{x — ttj), j = 1, . . . , K, and 9{x) the polar angle of x G M^. 

To push further the analysis, in [13] the attention is restricted to symmetric vortex config- 
urations in order to reduce the number of independent variables in E[a). In particular, the 
simplest ^—symmetric vortex configurations (a, n) (which are invariant under ^—rotations 
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and reflections w.r.t. the real axis) have the form: ao = 0, ai, . . . ,aN are the vertices of a 
regular A^— polygon with oi = 1 and rii = ■ ■ ■ = rij^f = m. We impose also the forceless 
condition Vi/(a) = 0, which simply reads as uq = —^^^m. 

Since |ai| = ■ ■ ■ = |aAr|, the only variable is the size a = |ai| of the polygon, and then the 
intervortex energy will be in the form E{a). Since A{a) is homogeneous of degree —2, we have 
that A(a) = ^, where 

2-Ki 2-Ki(N-l) 

Ao := A(l,eT^,...,e^^) (1.4) 

is the correlation coefficient for given no = —^Y^m and ni = ■ ■ ■ = un = m. In [13] the 
existence of c.p.'s of E{a) is shown for the cases {N,m) = (2,2) and {N,m) = (4,2) by 
comparing E{a) for a small and large, and using the positive sign of Aq (the correlation 
coefficient has value 8n and SOvr, respectively). It is also conjectured [13] that the correlation 
coefficient has values which are integer multiples of ^. With a long but tricky computation, 
in the next section we will disprove such a conjecture by showing 

Theorem 1.1. The correlation coefficient in always vanishes: Aq = 0, for all N > 2 

and m G Z. 



Beside the role of Aq in symmetry-breaking phenomena for fll.ip in M^, as already discussed, 
we will also explain its connection with the Ginzburg-Landau equation 

f -An = Pm(1 - |up) in VL ,^ . 

y u = g on dVt ^ ' ' 

on a bounded domain VL for strongly Type II superconductors k — )■ +00, where g : dVL — )■ 
is a smooth map. 

The energy functional for (11.51) 

has always a minimizer Uk in the space H = {u & H^{Q,C) : u = g on dQ}. When d = 
deg g 0, hy [2l[T5l|16] we know that on simply connected domains Uk has exactly \d\ simple 
zeroes ai, . . . , a\d\ for k large, where (ai, . . . , a\d\) is a critical point for a suitable "renormalized 
energy" W{ai, . . . , a\d\)- The symmetry-breaking phenomenon here takes place, driven by an 
external mechanism like the boundary condition that forces the confinement of vortices in 
some equilibrium configuration. A similar result does hold [2] on star-shaped domains for any 
solutions sequence Uk of (11. 5p . Near any vortex a^, the function u(| + a^) behaves like [/„.(x). 

Once the asymptotic behavior is well understood, a natural question concerns the construction 
of such solutions for any given c.p. (oi, . . . , ax) of W, and a positive answer has been given 
by a heat-flow approach [8l|9], by topological methods [1] and by perturbative methods [5|[T^ 
in case ni = ■ ■ ■ = uk = ±1- In [H], page 12, it is presented as an open problem to know 
whether or not there are solutions having vortices collapsing as ^ — >■ 00, the simplest situation 

2 

being problem (ll.Sp on the unit ball B with boundary value go = 

( -Au = k^u{l - |np) in B ,^ . 

\ u = go on OB. 

It is conjectured the existence of solutions to (II. 6p having a vortex of degree —1 at the origin 
ao = and three vortices of degrees -|-1 at the vertices laj, aj = e^^^~^^ for j = 1,2,3, of a 
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small (/ << 1) equilateral triangle centered at 0. This vortex configuration is ^—symmetric, 
forceless and has "renormalized energy" 

W^(0 = -67rln3-67rln(l-/6) + O(/9), / > 0. (1.7) 

In collaboration with J. Wei, we were working on this problem. Inspired by [3], we were aiming 
to use a reduction argument of Lyapunov-Schmidt type, starting from the approximating 
solutions Uk for (11. 6p given by 

3 

i=i 

with / —J- and Ik +oo, where the function (fk is an harmonic function so that Uk \dB= 9o- 
The interaction due to the collapsing of three vortices onto gives at main order a term {lk)~^ 
with the plus sign, i.e. for some Jo > there holds the energy expansion 

EkiUk) = 47i\nk + I + ^W{l) + Mlk)-^ + o{{lk)-^) 

= ATr\nk + I-3n\n3 + 3nf + Jo{lky^ + o{l^ + {lky^), (1.8) 

in view of (II. 7p . The aim is to construct a solution Uk in the form Uk[r]{l + ip) + {1 — r])e^], 
where ifj = ip{k) is a remainder term small in a weighted L°°{B)— norm and / = l{k) as 
k — >■ +00. The function is a smooth cut-off function with = 1 in U^^Qi?i/fc(/aj) and rj = 
in B \ U^^Qi?2/fc(^aj)- The function ip = ifji^) is found thanks to the solvability theory (up to 
a finite- dimensional kernel) of the linearized operator for (11.61) at Uk as / — ?■ and Ik — j- +oo, 
and by the Lyapunov-Schimdt reduction the existence of l{k) follows as a c.p. of 

Ek := Ek [UkW + + (1 - r,)e^^^^]) . 

If Uk is sufficiently good as an approximating solution of (II. 6p . we have that Ek = Ek{Uk) + 
o{{lk)~'^). Since Svr/^ + J^ilk)'"^ has always a minimum point of order k~i as k ^ +oo, by 
(II. 8p we get the existence of I = l{k) in view of the persistence of minimum points under small 
perturbations. 

Unfortunately, this is not the case. Pushing further the analysis, we were able to identify the 
leading term ipo = ipo{k) of i/j = ■0(/c), and compute its contribution into the energy expansion, 
yielding to a correction in the form: 

Ek = 4nlnk + I+ ^W{1) + Ji{lk)-^ + o{{lk)-^). (1.9) 

By (ll.7p and (II. 9p a c.p. l{k) of Ek always exists provided Ji > 0. First numerically, and then 
rigorously, we were disappointed to find that Ji = 0. 

Later on, we realized that — Ji is exactly the correlation coefficient Aq in (II. 4p (with = 3 
and m = 1) introduced by Ovchinnikov and Sigal [13]. If m is a solution of (ll.6p with vortices 
Co = and laj, aj = e^^^~^^ for j = 1,2,3, with uq = —1 and rti = n2 = = 1, then the 
function U{x) = ii(|) does solve 

-AU = U(l-\Uf) inB. 

U = go on dBk ^ ' 

with vortices oq and Ikaj of vorticities no = — 1, ni = ni = ns = 1. Since (11.11) and (ll.lOp 
formally coincide when k = +oo, it is natural to find a correlation term in the energy expansion 



SOME REMARKS CONCERNING SYMMETRY-BREAKING FOR THE GINZBURG-LANDAU EQUATION 5 

Ek in the form — ^ = Ji(/A;)~^, where a = Ik is the modulus of the Ika/s for j = 1, 2, 3. Even 
more and not surprisingly, the function Uk{f), where 

is a very good approximating solution for (11 .Gp which improves the approximation rate of Uk, 
does coincide with the refined test functions used by Ovchinnikov and Sigal [13] to get the 
improved upper bound (11 ■3p . 

In conclusion, the vanishing of the correlation coefficient Aq does not support any conjecture 
concerning symmetry-breaking phenomena for (II. ip or the existence of collapsing vortices for 
(II. 6p when k — ?■ +oo. Higher-order expansions would be needed in their study. 



2. The correlation coefficient 



Let > 2. Let 
for all J = 1, . . . , A^, ao = and tiq 



— ^ , J — 1, . . . , N, he the A^— roots of unity, and set rij = m E Z 
-^^^m. We aim to compute the correlation coefficient 
Aq = Ao(m) given in (11.41) . Since (in complex notation) V9{x) = \x\~'^{—X2, xi) has the same 
modulus as - = r^a, the correlation coefficient takes the form 



A, 



N 

E 

j=0 



N 

E 

j=0 



Hi 



(2.1) 



Since the integer m comes out as from the expression (12.11) . we have that Ao{m) = m'^Ao{l). 
Hereafter, we will assume m = 1 and simply denote Ao{l) as Aq. 

Let us first notice that Aq is not well-defined without further specifications, because the integral 
function in (12.11) is not integrable near the points aj, j = 0, . . . , N. Recall that the iV— roots 
of unity ai, . . . ,aN do satisfy the following sjTiimetry properties: 



N 



^4 = V|/|<A^, /^O, 



(2.2) 



as it can be easily deduced by the relation x^ 
is the validity of 



N 



Yl{x — ttj). A first application of (12. 2 p 



1 



N 

^x-a- 

j=i ^ j=i 



Nx 



N-l 



X 



N 



X^-1' 



(2.3) 



which implies that the integral function in (12. ip near has the form 
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N 



and is not integrable at when N = 2. Similarly, setting ak{x) 



2x 



k 



N, near we have that 



N 



j=0 ^ ~ 



N 

E 

i=o 



+ I 2Re — . T- h \ak[x)\ 

\x - afc r 



X - ak 
2 



-Re[(a: - afc)afc(a;)] 
16|x|4 



|afc(a:)|2 (2.5) 



a; — a.- 



The function can not be computed explicitly, but we know that 



2ai 



1 ^ 1 



-,N-l 



N 



1 " 1 

j=2 J 



(2.6) 



-.N-l 



N -1 



N 27rO-l) , . . 27r(i-l) 

1 — COS — -"^f^ + ^ sm — 



E 



N 



la 



N~ 



fl^2(l 



sm 



i=2 

27r0-l) 
TV 



COS 



27r0-l) - 



2(1 - cos 
= 



N ' 



in view of {a-,a^ ~ : j = 1, . . . , N, j ^ k} = {02, . . . , ctAr} and the symmetry of {ai, . . . , aj^} 
under reflections w.r.t. the real axis. By inserting (12. 6|) into (12. 5p we deduce that the integral 
function in (12. ip near ak has the form 



N 

T 

^ — ^ X — aj 



N 

E 

i=o 



X — ai 



\x - ak\ 



-Re[a'k{ak){x - at)^] + 0{- 



X - ak\ 



(2.7) 



and is not integrable at a^ when ^^(afc) 7^ 0. Since the (possible) singular term in (12. 4p . (12. 7p 
has vanishing integrals on circles, the meaning of ^0 is in terms of a principal value: 



An 



- lim 

4 e-S>0 



K2\uf^j,B,(afc) 



N 



E I 



N 

E 



rii 



(2.8) 



We would like to compute Aq in polar coordinates, even tough the set \ U^Qi?e(afc) is not 
radially symmetric. The key idea is to make the integral function in (12. 8 p integrable near any 
aj, j = 1, . . . ,N, by adding suitable singular terms, in such a way that the integral in (12. 8 p 
will have to be computed just on the radially symmetric set \ B^{aQ). To this aim, it is 
crucial to compute a^(afe). Arguing as before, we get that 



a'kiak) 



N-l 



N 



-E 



-.A'-2 



2 

2vr(j-l) N2 
N ) 



1 ^ 1 



^ f (1^ 
2 ^ 4(1- cos ^)2 

TV ^1 

^2(l-cos^) ' 



J=2 

2 27r(j-l) 

AT 



a,)2 



(2.9) 
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Af-1 TV 

Since there holds Yl ~ ^ ~ ~^ all A; = 2, . . . , in view of fl2.2p . we have that 
i=i i=2 



N 

n( 

i=2 



z'' - 1 



Af-1 N 



z — a,- 



p=0 j=2 



z — an 



p=0 



Z - flfc 



Ar-2 N-2-P 
p=0 «=o 



and then 



N N 

l[{l-aj) = N, l[{l-aj 

J=2 



J=2 



Ar-2 



^(iv-/-i)4. 



(2.10) 



By (1210|) we get that 



N 



AT 



4 ^ 4 



i=2 



A^2 
1=2 fc=2 



i=2 «,p=0 



"""'E''-4E'p-4E'--4(Eo-- ^'^-T'^-" ^(A--tf 



iV2 A^2 ^ JY2 



in view of Since by ([23]) ^'^(afc) = by (EH) we have that 



N N N 



j=0 ■> j=0 ■' 



^ (x — aj)2(l + |a; — ajP)' 



G L^(M'\{0}). 



Since 



lim 



lim 



we can re-write as 



Aq = - lim 

4 e^O 



K2\B,(0) 



(Ar + i)x^ + (iV_i) (AT -1)4 ^ 



2x(x^ - 1) 



16|x| 



X — ttj I 



|J^'[(x-a,)2(l + |x-a,f)- 



lim 



4 iR2\(B^(o)u{l-.<|x|<^}) 



16|x| 



-^Re 
4 



lim 

e-5-O 



(Ar + l)x^ + (iV-l),4 (AT-l)^ 
2x(x^ - 1) 



^ 1 



(o)u{i-e<H<^}) ~ + \x- Qjf ) 



ii-in (2.11) 



in view of (I2.3p . 
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As far as I, let us write the following Taylor expansions: for |x| < 1 there hold 
((A^ + l)x^ + (A^-l))2 



((AT _ 1)2 + 2{N^ - + (iV + 1) V^) ^(fc + 1)^'^^ 

fc>0 

{N -If + Y^ 'iNikN - = Ckx''^ (2.12) 

fc>0 



k>l 



and 



{{N - l)x^ + {N + 1)Y 



{1-x 



((A^ + 1)2 + 2{N^ - l)z^ + (A - 1) V^) ^(fc + 1)0;'=^ 

fc>0 

(A + l)2 + ^4iV(A;A + l)x'^^ = ^4a;'=^, (2.13) 



k>l 



k>0 



where = max{4A(A;A^ - 1), (A^- l)^} and 4 = max{4A(A;A^ + 1), (A+ l)^}. Letting e > 
small, by f l2.12p -( p.l3p we have that in polar coordinates (w.r.t. to the origin) I writes as 



I 



2-K 



pdp / d6 



16p4 



fc>0 



+ / pdp d9 



2n 



■kN^-ikNe\2 



^^-$:i$:(*+i)«r-ve*T 

(N-l) 



4 , N 



k>0 



■ike\2 



j=l k>0 



with Oe(l) — )■ as e — )■ 0, in view of 

= \af-'x-l\-^=\J2ik + l)a 



\x — a,- 



k{N-l) ki2 

j J- I , 



|1 — ajx\ 



^(fc + l)a; 



^ 1 2 
I 



A:>0 



fc>0 



for |x| < 1. By the Parseval's Theorem we get that 



I = 2n 



16 



l]|c.|V'=^-^-Aj](fc + l)V 



2k+l 



k>l 



k>0 



dp 



+2tt 



27rA 



1 

1-e 



k>l 



16p3 



A^(fcAr + A-1)V^' 



2kN+2N-3 



Y,{k+ifp 



fc>0 
2fc+l 



dp + o,{l) 



k>0 



k>0 



dp 



-2nN 



Ar^(A;A + Ar+l)V'''^"''^^' 



J](A; + 1)V''"' 



fc>0 
l-e 







+o,(l) = 27rA^ 
-2^(A; + l)V'=+i 



A^(A;A^ + A^- 1)V' 



fc>0 
2feAf+2iV-3 



TT 



4 + A^(A^' + l)- 



^^^(A;A^+A^+l)y'=™+l 



fc>0 



fc>0 



k>0 



TT 



dp + A(A' + l)- + o,(l) 
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as e — 7- 0. We compute now the integrals and let e — ?■ to end up with 

I = 2'kN ^^(A;iV + iV-l)p2'=™-2 + ^^(A;Ar + iV + l)p2fc™+2_^(^ + l)p2/c+2 

A:>0 fc>0 fc>0 

+N{N^ + l)'^. 

Denoting the function inside brackets as /(p), we need now to determine the explicit expression 
of /(p) for p < 1: 

fc>0 fc>0 fe>0 



2 ^ (l-p2^)2 2^ l-p2^ (1-p')' 

^2^2iV-2(l + p4) _ Arp2^-2(i _ ^4)^1 _ ^2W) _ 3^2^ ^ ^2i)2 

2 (l-p2^)2 

and then by the THopital's rule we get that 

N{N - l)p^-i + A^(A^ + l)p^+i - 2p( P'f + iVp'^"' - A^p'^+' 
4iVV(l) = 21im 

p-^i (1 — pY 

N-l 

-N'{N - 2)p^-2 - Ar2(iV + 2)p^ + 2( E /^T + M E />^')( E (j + ^)P^) 

j=0 j=0 j=0 

= hm 

P^i 1 — P 

(2iV + l)p^^ - (2iV - l)p^^-^ - p^-^ - p^ _ N\N' + 5) 
p->i 1 ~ P 3 

In conclusion, for I we get the value 

I = -Ar(Ar2 _ 1). (2.14) 
3 

Remark 2.1. In [13] the value of Aq was computed neglecting the term II in (12. lip . By (I2.14p 

notice that ^1= ^m^NlN"^ — 1) does coincide with Stt when {N,m) = (2,2) and SOtt when 
[NjUi) = (4,2), in agreement with the computations in |13j . 

As far as II, let us compute in polar coordinates the value of 

limV" / \ ^ = lim / pT{p)dp, 

where the function F is defined in the following way: 

N ,2. ^2 

^ de 



(pe*^ - aj)2(2 + p2 - ajpe-'(> - a^'^pe'^^ 
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with 7 the counterclockwise unit circle around the origin. Since 



w 



2 + p2 



2 + p 



,2 \ 2 



2p ' 

observe that w'^ — ^^j-ajw + a| vanishes at p±aj, with 

2 + p2 



,2 + p2 
• 2p 



P± 



2p 



± 



,2 + /92 



2p 



satisfying p_ < 1 < p+ in view of > V2. Since 



2 _ 



by the Cauchy's residue Theorem the function T{p) can now be computed explicitly as 



N 



dw 



P ,=1 ^7 

2nN 



f (PP- - 1) - ^ if p < 1 

\ {pp. - - pp-)-^ - p2 if p > 1. 



Since we have that 



(pp_ - 1)2 = ^(p2 - Vp4 + 4)2 = i(p4 + 2 - pVp' + 4), PP+ - PP_ = a/7+4, 



p^ + 2 p2 



we get that 

(pp_-l) (PP+-PP-) =-^^^^^^^^/___^^^^^^^^ 2V7T4 ' 2 

and the expression of r(p) now follows in the form 

p^ + 2 



r(p) = vrA^ 



v/7+4 



ttA^P^ 



Note that 



27rA^p2 if p < 1 
ifp>l. 



4p 



(2.15) 



v/7+4 

is integrable in (0, 00), and we have that 



p^ + 2 
v/7+4 



(p4 + 2)v/7T4 + p2(p4 + 4) 
^ s2 + 2 



1 r 

p^)dp = lim - / 



lim 

M-5-1-00 



M 



M2 



4 



Vs2 + 4 
M 



s)ds 



(2.16) 



lim — (VM2 + 4-M) = -. 
4 2 



Thanks to fl2.15p - fl2.16p we can compute 



lim 

e->0 



'(0,1-€)U{^,+«d) 

and for II we get the value 



/ pr(p)6?p = / pT{p)dp = 7rN, 

>/(0,l-e)U(T^,+oo) Jo 



II = -iV(A^2 - 1). 



(2.17) 



SOME REMARKS CONCERNING SYMMETRY-BREAKING FOR THE GINZBURG- LANDAU EQUATIONll 



Finally, inserting fl2.14p and fl2.17p into (12. lip we get that the correlation coefficient vanishes: 
Aq = 0. Then, there holds Ao(m) = for all m G Z, as claimed. 
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